Penetration depth study of superconducting gap structure of 2//-NbSe 2 
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We report measurements of the temperature dependence of both in-plane and out-of-plane penetration depths 
(A a and A c respectively) in 2//-NbSe2 . Measurements were made with a radio-frequency tunnel diode oscillator 
circuit at temperatures down to 100 mK. Analysis of the anisotropic superfluid density shows that a reduced 
energy gap is located on one or more of the quasi-two-dimensional Nb Fermi surface sheets rather than on the 
Se sheet, in contrast to some previous reports. This result suggests that the gap structure is not simply related 
to the weak electron-phonon coupling on the Se sheet and is therefore important for microscopic models of 
anisotropic superconductivity in this compound. 



The transition metal dichalcogenide 2//-NbSe2 is un- 
usual in that an incommensurate charge density wave state 
ffcDW < 33 K) coexists with superconductivity (T c ~ 7.1 K) 
til]. Recently, there has been renewed interest in this com- 
pound because its unusual superconducting properties, such 
as the field dependence of the thermal conductivity in the vor- 
tex state |0], bear a striking resemblance to the two-gap su- 
perconductor MgB-2- Many experimental probes have found 
evidence for strong gap anisotropy in 2H-NbSe2, but it re- 
mains unclear whether this results from there being different 
gaps on each Fermi surface sheet or some other form of gap 
anisotropy. Determining the nature of the gap anisotropy is 
an important step toward understanding its microscopic ori- 
gin. Is it related to the CDW, or does it rather result from an 
anisotropic electron-phonon interaction as in MgB2 OA? 

Angle resolved photoemission spectroscopy (ARPES) mea- 
surements yj] performed at 0.8 T c have found that the super- 
conducting gap is very small (below the experimental resolu- 
tion of 0.16 k-QT c ) on the small, pancake-like Fermi surface 
sheet associated with the Se p-orbitals. The gaps on the quasi- 
two-dimensional sheets, derived from the Nb d orbitals, were 
much larger (1.6 ± 0.16/cbT c ). The strength of the electron- 
phonon interaction (A ep ) is also different between these sheets 
Ojla], being rather weak on the Se sheet. This result suggests 
a similarity to MgB2 where a smaller superconducting gap on 
the 7r sheets [3] compared to that on the a sheets results from 
the being a large difference in A ep between these sheets com- 
bined with very weak a — ir scattering. The weak interband 
scattering is the crucial factor which produces the possibly 
unique, two-gap superconducting state in MgB2. It is as yet 
unclear whether a similar situation exists in 2//-NbSe2. 

We have performed measurements of the temperature de- 
pendence of the superfluid density to investigate the gap struc- 
ture in 2//-NbSe2. In an anisotropic superconductor, the 
temperature dependent components of the superfluid density 
Pi = l/\f(T) (i — a,c) are sensitive to the distribution of 
the size of the superconducting gaps on the Fermi surface |0]. 
Here we present measurements of the temperature dependence 
of both the in-plane and c-axis magnetic penetration depths 



[A Q (T) and A C (T)] of NbSe 2 from 100 mK up to T c . We find 
that our results are not compatible with the superconducting 
gap anisotropy as determined by the ARPES experiment [4]. 
Our data show that there is a reduced energy gap on one or 
more of the Nb Fermi surface sheets, but the gap on the Se 
sheet is not smaller than this, and may even be larger. 

In order to correctly interpret our measurements it is essen- 
tial to have a good understanding of the normal state electronic 
structure. Band structure calculations within the local density 
approximation J5[ S] show that three bands cross the Fermi 
level, giving rise to five Fermi surface sheets. Near the T point 
there is a small pancake-like sheet which derives mostly from 
the Se p bands (band 16) 0]. The other four sheets derive 
from the Nb d bands and are weakly warped tubes running 
along the c-axis, centered on the T and K points. The two 
surfaces derived from the bonding Nb d band (17a and 17b) 
are significantly more warped than those from the antibonding 
Nbd band (18a and 18b). 

This band structure has been found to be in good over- 
all agreement with ARPES J3, 0] an d de Haas-van Alphen 
(dHvA) measurements [5]. Significantly, the formation of the 
CDW state does not seem to lead to major Fermi surface re- 
construction. dHvA measurements have only resolved signals 
originating from the small Se pancake and show that, although 
the shape is similar to bandstructure calculations, its size is 
somewhat smaller. The dHvA and ARPES results show that 
the mass renormalization factor (1 + A m * ) varies considerably 
between sheets. It is largest on the Nb sheet 17b (A m » ~ 1.9 
from ARPES §6]) and smallest on the Se sheet 16 (A m * ~ 0.3 
from dHvA H]X On the other Nb sheets A m « ~ 0.85 0]. 

In the clean local limit, pi (0) is related to the unrenormal- 
ized plasma frequencies uj p .i (in eV) by p s — (u> p ^e/ch) (1 + 
A m *) _1 CD- Band structure calculations of oj Pt i show that 
the Se sheet contributes only ~ 2% to p a (0), but ~ 80% to 
p c (0) (this reduces to ^50% if the Se band is shifted slightly 
to agree with the dHvA results (§]). This shows that p a {T) is 
only sensitive to the gap on the quasi-2D Nb cyclinders, while 
p c (T) is sensitive to the gap on both the Se pancake sheet and 
the Nb sheets in comparable proportions. 
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FIG. 1: Temperature dependence of the in-plane penetration depth, 
AA a in single crystals of NbSe2. Data for three different samples are 
shown, offset for clarity. The absolute scale of AA a is most accurate 
for sample 1 (see text). For sample 2 the data have been divided by 
1.3. Inset: Frequency shifts at higher temperature showing the charge 
density wave transition. 



The temperature dependence of the penetration depth was 
measured in single crystals of NbSe2 using a radio frequency 
(RF) resonant LC circuit driven by a tunnel diode [9], oper- 
ating at ^12 MHz and low magnetic field (-£/rf — 10~ 2 Oe 
<C H c i). The sample was attached with vacuum grease to 
a sapphire rod, the other end of which is attached to a tem- 
perature controlled stage. The sample is placed in a copper 
solenoid which forms the inductor in the LC circuit. Changes 
in the resonant frequency of the circuit (AF) are directly pro- 
portional to changes in the field penetration in the sample. 

For H ||c the screening currents flow only in the basal plane 
and AF = aA\ a , where the geometric factor a is estimated 
using the technique described in Ref. II 1 CHI - For H±c currents 
flow both in the plane and along the c-axis, and (neglecting 
the small contribution from the top and bottom faces) for a 
rectangular sample with dimensions l x , l y and l z [x,y are in- 
plane and l x ,l y ,l z 3> A] AF for H \\y is given by itlOll : 



AF/AFq = 2AX a /l z + 2A\ c /l x . 



(1) 



Here AFo is the frequency shift obtained when the sample is 
completely removed from the coil and accounts for the sample 
demagnetizing factor as well as the coil calibration factor. 

Measurements were conducted on samples from three dif- 
ferent sources (Tsukuba, Lausanne, Bell labs) in three differ- 
ent laboratories. All samples were grown via the usual iodine 
vapor transport technique and are known to be of high qual- 
ity, having a high T c ~ 7. IK and low residual resistances 
(RRR~ 40-80). Some of the samples were from the same 
batch as those used for dHvA measurements Q5D- Experiments 
in Bristol were performed in either a dilution fridge (F m i n ~ 
100 mK) or a pumped 4 He cryostat (T m [ n ~ 1 .3 K). In Greno- 
ble and Urbana-Champaign a 3 He cryostat was used (T m ; n ~ 
0.5 K). In total more than 15 samples were measured. 

Fig- [U shows the low temperature variation of the in-plane 
penetration depth, which, as mentioned above, probes mainly 



the excitations on the Nb sheets. Data are shown for three 
samples measured in different laboratories. Sample 1 (Bristol) 
was very thin (0.48 x 0.55 x 0.011 mm 3 ) with HLc, hence 
the contribution from AA C is negligible (< 2%). Samples 2 
(Urbana/Ames) and 3 (Grenoble) were thicker, having l x ~ 
l y = 2.3l z , and were measured with 7F||c. All three curves 
are fitted to the expression 
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which approximately reproduces the full solution to the BCS 
equations for T < F c /3. For all samples we find Ao = 
1.1 ± 0.1 kftT c (the quoted error includes a weak dependence 
on fitting range). In a multigap superconductor, or one with 
a distribution of gap values, at low temperature Ao is a mea- 
sure of the smallest gap. For example, in MgB2 a low tem- 
perature fit measures the gap on the it bands, which is much 
lower than the BCS weak coupling value 111 111 . Combined 
with our knowledge of the normal-state electronic structure, 
the data clearly show the presence of excitations with an en- 
ergy gap much smaller than the weak coupling BCS value 
(Ao = 1.76 k^,T c ) on the quasi-2D Nb sheets of Fermi sur- 
face. 

The temperature dependence of AA a was very similar in 
all samples. We find that the absolute values of AA a mea- 
sured with H\\c are sometimes higher than those for HLc 
(e.g., sample 2 in Fig. Q]). This is in the opposite sense to 
that expected from the additional A c contribution for HLc. 
It is likely that this results from the mica-like morphology of 
the crystals. Although very flat (001) faces may be prepared 
by cleaving, cutting the crystal perpendicular to this direction 
produces splintered edges which may have a larger effective 
area than their geometric cross-section (and hence larger ef- 
fective field penetration). For this reason, we believe the cali- 
bration factor is most accurately determined for the thin sam- 
ple with HLc (sample 1 in Fig. [TJ where ^98% of the sig- 
nal comes from the flat faces and demagnetization effects are 
small. 

We now turn to an analysis of A c which probes the gaps 
on both the Se and Nb sheets. If the gap on the Se sheet is 
significantly smaller than that on the Nb sheets J4J we would 
expect the temperature dependence of A c to differ markedly 
from that of A a . 

For thick samples with HLc a significant proportion of the 
total measured AF(T) comes from currents running along the 
c-axis. In the inset to Fig. [2] we show data for sample 2 where 
l z /l x ~ 0.44. AF is temperature independent for T < 1.2 K. 
This clearly shows that the gap is not substantially smaller on 
the Se sheet than on the Nb sheets. 

By measuring AF with H||ir before and after cutting the 
sample in half along the field direction (see Fig. [2]) we are 
able to make a more precise determination of AA C (T). In 
principle, the in-plane contribution from the large faces is un- 
changed after cutting, so the difference between frequency 
shifts is due to c-axis screening currents on the newly exposed 
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FIG. 2: Low temperature behavior of AA C determined by cutting the 
sample (see text). Data for two different samples are shown (samples 
4 and 5). For sample 5 the data have been divided by 2 and offset for 
clarity. The results for AA a (T) (sample 1) are shown for comparison 
(dashed line). The solid line show a fit to Eq. l[2} up to 2.4 K. The 
lower inset shows the geometry of the samples before (I) and after 
(II) cutting. Upper Inset: Raw measured frequency shift for a thick 
sample with HJlc (same as sample 2 in Fig.[T). 



surfaces [Eq. ([T)]. In practice, there may remain some contri- 
bution from AA a because the effective edge area may change 
when the sample is cut due to splintering. However, we found 
that when the sample was cut in two stages (first in half then 
in quarters along the field direction), each cut resulted in an 
identical additional contribution to AF(T). 

In the main part of Fig. [2] we show our result for the 
low temperature behavior of AA C determined by the cutting 
method in two samples AA C has a very similar temper- 
ature dependence to AA a . There is a variation in the absolute 
magnitude of AA C between samples, possibly due to a stray 
contribution from A a , but the measured temperature depen- 
dence of A c is very similar. In this temperature range AA C (T) 
displays exponential behavior at low temperature indicating a 
minimum gap, Ao = 1.2 ± 0.2 T c , similar to that found in 
AA (T). The absence of any strong temperature dependence 
in AA C (T) at low temperature and the similarity to AA (T) 
demonstrates our key result, that the gap on the Se sheet is at 
least as large as that on the Nb sheets. 

In Fig. [3] we show the normalized in-plane superfluid den- 
sity, p a = A a (0) 2 /A a (T) 2 for a range of values of A a (0) close 
to that determined by muon spin rotation measurements lfl3ll 
[Aa(0) = 1250 A]. As the data clearly do not follow the stan- 
dard BCS behavior we have attempted to model p a (T) using 
a fc-dependent gap. In the clean local limit BCS theory p a is 
given by 07D 
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where dSp is an element of Fermi surface, v is the Fermi 
velocity, / is the Fermi function and A^ is the k dependent 
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FIG. 3: Temperature dependence of the normalized in-plane super- 
fluid density p a calculated using the A a (0) values indicated. The 
solid lines are fits to the two gap model described in the text. 
The dashed line is the weak coupling BCS prediction with Ao = 
1.76 k B T c . 



energy gap. 

Our first approach is to assume that the gap on one, or more, 
of the two quasi-2D Nb sheets is different to the others, but 
there is no variation of the gap within each sheet (similar to 



the situation in MgB 2 ) till [MO]. Eq. © then reduces to 



p(T) = xp[T, Ai(T)] + (1 - x)p[T, A 2 (T)] 



(4) 



where x is the Fermi surface weight, which is proportional to 
the total superfluid density on sheet(s) with the reduced gap, 
and Ai, 2 (T) = A 1 , 2 A B cs(T)/A B cs(0), where A BC s(T) is 
the BCS weak coupling gap. The fits to p a (T) with this model 
are excellent (see Fig. [3). The parameters are given in TableQ] 
The value of the small gap varies little with A a (0) and is close 
to that found from fitting Eq. (O at low temperature. The gap 
ratio is A 2 /Ai = 1.8 ±0.2. x ~ 0.5 indicating that the sheets 
with the reduced gap contribute ~ 50% to the total in-plane 
superfluid density. 

Our second approach considers a gap with 6-fold symmetry 



TABLE I: Variation of the fit parameters to models of the in-plane su- 
perfluid density as a function of the assumed zero temperature value 
of A a . Ai, A2, and x are the parameters of the two-gap model [Eq. 
10] and A(0) and e are the parameters of the 6-fold anisotropic gap 
model [Eq. J5)]. The final line shows the parameters found by fitting 
the same models [14] to the heat capacity C data of Ref.ll5L 



A„(0) # % 



6-fold gap 

A mln (0) 



l+e 
1-e 



Two gap 

Aa 

T c Ai T c 

1200A 0.99 1.59 0.43 0.94 1.74 

1400A 1.04 1.76 0.49 0.91 2.07 

1600A 1.08 1.91 0.50 0.91 2.33 



C 



1.31 1.76 0.30 1.25 2.31 
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(NbSe2 is hexagonal) on all of the Nb sheets, 

A(^T) = A min (T) [1 + ; COs(6 ^ (5) 
1 — e 

where, as above, A m ; n (T) is assumed to have the BCS tem- 
perature dependence and <f> is the in-plane angle (we approxi- 
mate the Nb sheets as simple tubes). We find that this model 
fits PaiT) as well as the two gap model above. The two free 
parameters [A min (0) and e], which vary slightly with the as- 
sumed A a (0), are given in TableQ] The ratio of the maximum 
to minimum gap is (l + e)/(l — e) = 2.0 ±0.3, which is close 
to that found for the 2 gap fits. 

A two-gap analysis of p c {T) has also been performed. The 
minimum gap is found to be Ai ~ 1.2 ± 0.2 fcsT c , in agree- 
ment with the low temperature fit of AA C (T) to Eq. ©. The 
value of the larger gap ranges from 1.2 to 1.8 fc^Tc depending 
on our assumed value of A c (0) and the relative contribution of 
the Nb and Se sheets to p c (T). 

By analyzing the zero field specific heat data of Ref. 
with the same models iflrjll we find that the gap values are 
~ 30% higher than those from the superfiuid density analysis. 
However, the values of the gap anisotropy are very similar (see 
TableHJl lH. 

The above analysis is unable to differentiate between a dis- 
crete two-gap model and a continuously varying gap. How- 
ever, the influence of impurities strongly favors the two-gap 
model. We find no difference in the behavior of A a for sam- 
ples with RRR in the range 40-80. A similar conclusion was 
found from analysis of the field dependence of the specific 
heat 111 (ST . The criterion for observing an anisotropic (or sheet 
dependent) gap is Ht^ 1 <C \J (A)SA, where t _1 is the im- 
purity scattering rate, (A) is the average gap and SA is its 
variation over the Fermi surface [1911 . Using the parame- 
ters of either our fits yields \J {A)5A ~ 7K. For a sample 
with RRR of 40, we estimate the isotropic mean free path is 
~ 270 A, an d for -1 = hv a /(k B £) ~ 27 K which is ~4 times 
larger than y/ (A)5A. The scattering must therefore be highly 
anisotropic. A likely scenario is that there is reduced scatter- 
ing between the two Nb bands. This would allow two differ- 
ent gaps to exist on these pairs of Fermi surface sheets (as for 
MgB2 [20]). The reduced energy gap on one of the Nb sheets 
could be related to the CDW state formation. 

In conclusion, our measurements of the anisotropic super- 
fluid density in NbSe2 show evidence for anisotropic, proba- 
bly multi-band superconductivity. However, this is not of the 
form suggested by some previous reports, where the small Se 
sheet is the only sheet with a reduced gap, and is therefore not 
simply due to the weak electron-phonon coupling on the Se 
sheet. We hope our results will stimulate theoretical investi- 
gations of the microscopic origin of the gap anisotropy, which 
may be related to the formation of the CDW state. 

We thank I. Mazin for useful discussions, as well as H. 
Berger, Y. Onuki and R Gammel for providing the NbSe2 
samples. 
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